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Obtaining cross sections for nuclear reactions at intermediate energies based on the Glauber 
formulation has a long tradition. Only recently the energy regime of a few hundred MeV has 
become accessible to ab-initio Faddeev calculations of three-body scattering. In order to go to higher 
energies, the Faddeev equation for three-body scattering is formulated and directly solved without 
employing a partial wave decomposition. In the simplest form the Faddeev equation for interacting 
scalar particles is a three-dimensional integral equation in five variables, from which the total cross 
section, the cross sections for elastic scattering and breakup reactions, as well as differential cross 
sections are obtained. The same observables are calculated based on the Glauber formulation. The 
first order Glauber calculation and the Glauber rescattering corrections are compared in detail with 
the corresponding terms of the Faddeev multiple scattering series for projectile energies between 
100 MeV and 2 GeV. 

PACS numbers: 21.45-v,24.10.Ht,25.10.+s 



I. INTRODUCTION 

The formulation of high energy scattering introduced by Glauber [l|, [^j has a long tradition in being applied to 
calculating nuclear reactions with hadronic as well as electromagnetic probes. Today Glauber amplitudes are widely 
used to account for final state interactions in (e,e'p) reactions at high and intermediate energies [1, 0, [1]. Glauber 
theory is applied to investigate color transparency, in {e,e'p) reactions and in the electroproduction of mesons @, 01- 
In addition, it is applied to heavy ion reactions [8|. Most recently, several generalizations of Glauber theory have 
been employed for the description of electron scattering mainly from lighter nuclei, e.g. the generalized eikonal 
approximation (GEA) [1, [l^, or relativistic eikonal approaches [Ullilj. These methods are applied both to A{e,e'p) 
and A{p, 2p) reactions. Nuclear scattering at intermediate energies in the context of radioactive beams also takes 
advantage of this formulation [3, [ij] . 

Despite the widespread use of the Glauber formulation in reactions involving few as well as many-body systems 
there has been little work on rigorous tests of the accuracy and/or limits of the Glauber ansatz, apart from some 
studies concerning the quality of the eikonal approximation for (e, e'p) reactions as well as for elastic proton 
scattering from halo nuclei [Tgl . 

Only recently exact Faddeev calculations for three-body scattering in the intermediate energy regime became 
available. This progress is based on a formulation of the Faddeev equations, which is based directly on momentum 
variables and does not rely on traditional partial wave expansions. The formulation and numerical realization for the 
nonrelativistic Faddeev equations [l3, ^ well as fully Poincare invariant ones [l^, [10, HH have been carried out 
for scalar interactions up to projectile energies of 2 GeV. This now allows us to perform a detailed comparison of 
three-body Faddeev calculations of total and elastic cross sections with the ones obtained by Glauber calculations, 
and thus test the range of validity of Glauber calculations at least in the three-body system. 

The naive expectation is that for sufficiently high energies the total n+d cross section is just the sum of the two total 
cross sections for two-nucleon scattering. The question whether this is valid is an old one and has been investigated 
in the context of the eikonal approximation [l[ [2^, [l^, . It turns out that rescattering (shadowing) corrections 
need to be added. Since at the time no exact calculations were available, the question remained open, whether those 
corrections were sufficient and at which energy the eikonal approximation became valid. We will investigate such 
questions in this work and mostly follow the formulation given in Refs. [25l . [26j . In order to make closer contact 
with that work as well as for transparency, we will use the nonrelativistic formulation of the Faddeev equation for 
our comparison, even though we carry out calculations in the intermediate energy regime. We will also use the 
same two-body scattering amplitude as well as deuteron wave function in the Glauber calculation as enters in our 
Faddeev results, so that we can clearly identify the effects of the Glauber approximation on the multiple scattering. 
In Section II we will briefly present the Faddeev framework and its multiple scattering expansion, and in Section III 
we will re-derive the essential expressions of the Glauber approximation necessary for the comparison to the Faddeev 
formulation. Our numerical results for the total and differential cross section for elastic scattering in both formulations 
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and their discussion will be given in Section IV. We summarize and conclude in Section V. 



II. THE FADDEEV MULTIPLE SCATTERING SERIES 

Various presentations of three-body scattering in the Faddeev scheme are presented in the literature [22|, [2^ . We 
solve the Faddeev equation for three identical particles in the form 

T\<j>)^tP\(l>)+tPGoT\(j>). (2.1) 

The driving term of this integral equation consists of the two-body t-matrix t, the sum P = P12P23 + ^'13^23 of a 
cyclic and anti-cyclic permutation of three particles, and the initial state 10) = I^^Qo), composed of a two-body bound 
state (j>d and the momentum eigenstate of the projectile particle. The kernel of Eq. (|2.1|) contains the free three-body 
propagator, Gq = {E — Hq + ie)~^, where E is the total energy in the center-of-momentum (cm.) frame. 
The operator T determines both, the full break-up amplitude 

C/o = (1 + P)T, (2.2) 

and the amplitude for elastic scattering 

U = PG„^ + PT. (2.3) 

In this paper we focus on three identical bosons and use a momentum space representation. For solving Eq. (|2.ip . 
we introduce the standard Jacobi momenta p for the relative momentum in the subsystem, and q for the relative 
momentum of the spectator to the subsystem. The momentum states are normalized according to (p'q'|pq) = 
6^{p' — p)i5'^(q' — q). Projecting Eq. (|2.ip onto Jacobi momenta leads to [l^ 

(pq|T|</)dqo) = (/-d ^q+ ^qo^ ^p, ^q + qo,-E- ^0^2 

V 2 Am J E - -^{q^ + q"^ + q- q") + i£ 

Here ts{p',p) ~ t{p,p') +t{—p',p) is the symmetrized two-body t matrix and the total energy E is explicitly given 
as 

E^Ed + -^ql =Ed + ^Eiab. (2.5) 

We assume that the underlying two-body force generates a t-matrix as solution of a two-body Lippmann-Schwinger 
equation, and that the force supports one bound state with energy Ed- Then Eq. (|2.4p is solved directly in momentum 



space as function of the vector Jacobi momenta and angles between them. For further details we refer to Ref. [17 1. 

In this work we want to study the high energy behavior of the cross sections. In the framework of the Faddeev 
formulation, we consider the multiple scattering series defined by the integral equation given in Eq. (j2.ip . For the 
elastic scattering amplitude this leads to the series expansion 

(*|J7|$) = ($|PGo-i|$) + ($|PtP|$) + ($|PtPGotP|$) + • • • , (2.6) 

which for the fully converged solution is summed as Pade or Neumann series. Via the optical theorem the total cross 
section is related to the imaginary part of the operator U as 

afS = -{27:f^<^m{^U\^). (2.7) 
3go 

This again can be expanded in a multiple scattering series through 

2i3m($|C/|$) = ($|J7|$) - ($|C/|$)* 

= ($|P(i- it)P|$) + ($|PiPGotP|$) - ($|P<^GSPttp|$) • •• (2.8) 



It turns out that the Faddeev multiple scattering series for the total cross sections converges fast for projectile energies 
higher than 500 MeV [13, H!!- This finding is encouraging for our study of the Glauber formulation. 
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One can also start from Eg . (12.81) and derive a high-energy hmit in the frame work of the Faddeev equations. This 
has been carried out in Ref. [2J|. Here the terms of first and second order in t have been investigated in the hmit 
E ^ oo leading to the analytic result 




= 2afJ^ + Oif)+Oit') (2.9) 

In first order this expansion gives twice the two-nucleon cross section. The second order correction terms can now be 
tested against exact Faddeev calculations. In Ref. [l^l this was not possible. 



III. THE GLAUBER CALCULATION 



For the convenience of the reader we briefly sketch the derivation of the Glauber amplitude for elastic n + d 
scattering, following the work given in Refs. [25|, [2^. In the laboratory system with the projectile momentum k and 
the momentum transfer q = k — k' the two- particle scattering amplitude according to Glauber [ij is assumed to have 
the form 



where r(b) is given through the eikonal phase x(b) as 

r(b) = 1 -e^xC'^ 



(3.1) 



(3.2) 



The vector b is assumed to be perpendicular to a direction n, specified below. Scattering off a deuteron target then 
leads to 



ik 
2^ 



dH e^i-'^ J d\ rd{r){l - e''(x(b-i/2s)-hx(b+i/2s))) ^^^^^^ 



(3.3) 



where (/)d(r) is the deuteron wave function and r = s -I- h{h ■ r). Using these definitions it follows that 



r(b - U) + r(b + is) - r(b - is)r(b + 



0d(r) 



Using the inverse of Eq. (|3.ip . 



ge-''i-^A.(q), 



(3.4) 



(3.5) 



and again integrating over a vector q in a plane perpendicular to fi one can eliminate the F's in favor of the two-particle 
scattering amplitude leading to the expression given in Refs. [2^, [2^ 



i^fc(q) = 25(iq)/(q) + ^ j dV^(q') A(^q + q')/fe(^q - q')- 



The deuteron wave function occurs in the form factor 



5(q) = j dV|0,(r)|2e^q- = j dV|0,(r)pe''i" 



From the above expression one obtains the laboratory differential cross section for n + d scattering as 



d^i 



\FM)Y 



(3.6) 



(3.7) 



(3.8) 
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and the total cross section for n + d scattering 



ND 



4tt 

— 3mFfe(0) 

fc 



(3.9) 



In first order in the two-body scattering amplitude Eq. p.9p leads to the expected result for the total cross section 

= 25(0)^3mA.(0) = 2^3mA.(0) = 2afJ^ (3.10) 
and to the differential cross section for clastic n + d scattering 



^^ND,l.t 

dn, 



-q IA(q)l 



For the total cross section for n + d scattering, including the second order correction, one obtains 



2a. 



NN 



Sat, 



where 



Sato 



^ J d'q' 5(q')/fe(q')/fc(-q')- 



(3.11) 
(3.12) 
(3.13) 



This correction term represents rescattering events. 

Adding the rescattering corrections to the differential cross section for elastic scattering leads to 



da 



ND 



dfli 



1 



—S ( -q ) 3m 



2 

irk 
1 

27rfc 



l/fe(q)P 

/fe(q) I d\' S{q') /fc(l/2q + q') /fc(l/2q- q') 
J dV ^(q') A.(l/2q + q') /fe(l/2q-q') 



(3.14) 



Next, we need to discuss the choice of the two-particle scattering amplitude /fc(q). In the original work of Refs. [l|, 
[Hjllll, the cikonal phase is given in terms of the potential V{r) as 



X(b) 



1 r°° 

- V(h + z) dz 

V J-oo 



V 

2m 



V{Vb^ + z^) dz, 



(3.15) 



Where z is a vector parallel to k. The above definition would result in a two-body scattering amplitude in the eikonal 
approximation. Studies of the quality of the eikonal approximation for potential scattering have been carried out 
in detail in the past for different potentials. The for us interesting case of Yukawa-type potentials is considered e.g. 
in Refs. [l^, [1^. Since our main interest is the comparison of the n + d cross sections, we need to use the same 
two-particle scattering amplitude /fc(q) for our Glauber calculation as we use for the Faddeev calculation. Thus we 
need to use the scattering amplitude obtained from the solution of a Lippmann-Schwingcr equation with the potential 
V as driving term. 

This two-particle scattering amplitude /c(p',p) is given in the cm. frame and is related to the on-shcU two-body 
t-matrix as 



(p'|y|^(+)) = {p'\t{E, = p!|p) ^t(p',p;5), 

and leads to the scattering amplitude 



m' 



/c(p',p) = -?(2^)2t(p',p;^) 



(3.16) 



(3.17) 



For a scalar potential i is also a scalar and depends on the magnitude of |p| = |p'| and the angle Xc =^ p' • p between 
the two vectors p' and p [29j . 
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i(p',p; — ) = tip^Py^c] — ) = t{j),xc). 

m m 



(3.18) 
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For identical particles we use as in the Faddeev calculations the symmetrized t-matrix ts{p',p; — ). The cm. differ- 
ential cross section for the scattering of two identical particles in then given as 



and the laboratory differential cross section as 



^ = l/cb,^c)P (3.19) 



da da , , 



where the standard relation between cm. and laboratory scattering angles for two equal mass particles is given as 



X/ = k • k' = ' - ' 

Xc = 2xf - 1. (3.21) 

If one defines 

^ = l/r(fc,x/)P (3.22) 

the two-particle scattering amplitude in the laboratory frame reads 

fnk.xi) = 2^if^{p,x,{xi)). (3.23) 

All that remains now is choosing the unit vector n. The original suggestion by Glauber [l|, [2^ is 

(a) : n = k, (3.24) 

with the approximation |k| = |k'|. In other works [25l. [27l. [28j the choice 

(b) : n = k + k' (3.25) 

is preferred. We will use both choices in our calculations to explore the differences. The momentum transfer q enters 
the scattering amplitude in the form of 

xe = l-^ (3.26) 

The Glauber ansatz requires to replace the momentum transfer by q^^ with q^ being orthogonal to n. For the 
two cases of Eqs. p.24p and (|3.25p this leads to 

(a) : q^2 ^ k\l - xf) 

(b) : q^2 ^ fc2(i„^2)_i^^ (3 27) 



For forward angles both choices become equivalent. 

The total cross section for two-body scattering entering into Eq. p.lOp is then explicitly given as 



a^J^ = J dn,\f!ip,x,)\'. (3.28) 

The scattering amplitude /fc(q) entering Eq. ()3.1ip can now be directly obtained from Eq. (|3.23p as 

fkici)^ ftik.xi). (3.29) 

In calculating the angle Xc{xi), the momentum transfer of Eq. (j3.26p must then be replaced by g^, and xi is related 
to Xc through Eq. (p?2T|) . 

For the second order correction to the differential cross section the integral in Eq. (|3.14p needs to be evaluated. 
According to the Glauber ansatz the argument of the scattering amplitude needs to be evaluated as l/2q_L ± q', so 
that the integration is carried out in the plane defined by n. Explicitly, the integral is calculated as 



/rfV 5('z')A-(^q + q')/fc(^q-q') 
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i ( {2n)^ ['"^^^ dq'q' S{q') r d<P ^ xi{xM)xi{x(~)) U{p, x^+^)ts{p, x^-''), (3.30) 



with the definition of ts{p,x) from Eq. (|3.18p and 

xW = l-iM±i4£^ii±^. (3.31) 

Here </> is the angle between the vectors qx and q'. The upper Umit Qmax is the maximum momentum transfer aUowed 
for a given projectile laboratory energy, corresponding to Xc = — 1 in the angle argument of the <-matrix. 
For the second order correction to the total cross section, Eq. (|3.12p . we only need to consider the special case q± ^ 0. 
Then the integral of Eq. (|3.30p simplifies to 



Satot = ^rn^27rf / ""^^ dq' q'S{q ')xi{x) 3?e [tl{p, x)] , (3.32) 



with X = 1- 



In Ref. [26| a further approximation is suggested arguing that if the form factor of the deuteron decreases much 
more rapidly than the scattering amplitude as function of q' , then the integral in Eq. (|3.32p may be approximated by 

Satot = ^m2(27r)5 SRe [tlip, x = 1)] /'""^ dq' q'S{q'), (3.33) 

Jo 

where the scattering amplitude in forward direction {x — 1 or 9 — 0) is taken out of the integral. 



IV. RESULTS AND DISCUSSION 



Our explicit calculations are based on an interaction chosen as superposition of two Yukawa interactions of the 
Malfliet-Tjon ^ type, 

in which the parameters, given in Table [H are fitted such that a bound state, the 'deuteron' is supported with a 
binding energy = —2.23 MeV. This interaction enters the nonrelativistic Faddeev equation for identical bosons via 
the symmetrized two-body i-matrix. The Faddeev equation is exactly solved without a partial wave decomposition, 
using momentum vectors and angles between them. This allows us to calculate three-body observables in the GeV 
region. The details of the calculations are given in Ref. The deuteron wave function entering the form factor of 
Eq. (|3.7p is obtained as solution of a bound state equation with the potential of Eq. (|4.ip as input. As a note, we are 
of course aware that this is a model study, in which the two-body i-operators are generated through a non-relativistic 
Lippmann Schwinger equation from a simple Hermitian two-body force. 



A. Total Cross Sections 



First we want to explore the convergence of the Faddeev multiple scattering series when considering total cross 
sections. The top panel of Fig. [1] shows the total cross section for three-body scattering calculated using Eq. (|2.7p as 
function of the order in the two-body i-matrix when adding up successively the terms of the multiple scattering series 
created by Eq. (j2.ip . The convergence is shown for selected projectile laboratory energies ranging from 200 MeV 
to 2 GeV. For energies larger than 1 GeV the curves are essentially fiat, meaning that the first order term given 
by T ^ tP is already sufficient to capture the result of a full Faddeev calculation. At 500 MeV one needs at least 
one iteration (or one rescattering term) to reach the exact Faddeev result, whereas at 200 MeV one needs more 
terms, which however give small contributions. The middle part of the panel shows the total cross section for elastic 
scattering and the bottom one the total cross section for breakup reactions. While the total cross section for elastic 
scattering converges very fast for projectile energies above 500 MeV, the one for breakup reactions requires at least 
one rescattering term even in the GeV range. 

Our Glauber calculation contains only a first order term and a rescattering correction. First, we concentrate on 
the total cross section and compare each order separately with the corresponding order of the Faddeev calculation. 
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The first order terms are compared in Fig. [2] as function of the projectile laboratory energy. According to Eq. (|3.10p 
the first order contribution to the total n + d cross section is simply given by twice the two-body total cross section. 
From Fig. [2] we see that for energies greater than 200 MeV both first order terms agree. This is not surprising, since 
the first order Faddeev term, T = tP, contains twice the two-body t-matrix, the permutation operator guarantees 
scattering contributions from both constituents of the target. Our Glauber calculations start from the same two-body 
t-matrix as the Faddeev calculation. From this point of view it is not surprising that both first order calculations 
agree so well. The same conclusion was already reached in Ref. [23 |. 

A more crucial test is a comparison of the second order correction term. In Fig. [3] we show the first rescattering 
correction to the n + d total cross section, which corresponds to second order contribution in t of the Faddeev 
multiple scattering series. The Faddeev result is shown as solid line. A first observation is that this rescattering term 
contributes significantly at lower energies, and is still present at 2 GeV. Thus, the very naive expectation that the 
n + d total cross section comprises only the scattering from the target constituents is not fulfilled in the energy regime 
up to 2 GeV which is considered here. There is always a contribution due to rescattering. The second order Glauber 
correction Sa from Eqs. (|3.13p and (|3.30p is given by the short dashed line. Here we see that in contrast to the first 
order contribution the two lines only start to agree around 1 GeV. This may be understood when having in mind that 
the second order Faddeev term is explicitly given as tPGotP, where the free propagator Go describes the arbitrary 
free motion of the three particles in the intermediate state. The Glauber second order term assumes that the target 
constituents are frozen. At sufficiently high energy this assumptions seems to be sufficient to describe the rescattering 
within the target, which sometimes is referred to as 'shadowing' correction. 

Franco and Glauber [2^ suggest a further approximation to Sa. They argue that if the deuteron form factor 
decreases much more rapidly than the scattering amplitude as function of the integration variable q' , then the integral 
of Eq. p.32p may be approximated by the expression of Eq. (j3.33p , where the scattering amplitude is approximated 
by its value in forward direction and taken out of the integral. The double-dash-dotted curve of Fig. [3] represents this 
approximation. It is quite obvious that this approximation is too simplistic and gives a second order contribution 
almost double the size of the original correction at all energies considered. The long-dashed curve represents the 
analytic evaluation of the second order Faddeev term in the optical theorem of Eq. (12.91) . The contribution calculated 
in there overestimates the second order contribution in a similar fashion. In Ref. [2J|, integrals over the two-body 
t-matrix and the deuteron wave function were evaluated using the method of steepest descent. This is equivalent to 
extracting the dominant part of an integral as a constant from the integral. The failure of those two approximations 
should lead to the conclusion, that although it is true that for higher energies the scattering amplitude is peaked in 
forward direction [2^, the integration over the form factor is important and should not be approximated further. 

After comparing the terms of the multiple scattering series separately, we show in Fig.|4]the sum of the contributions 
as function of the projectile laboratory energy. The solid line represents the exact Faddeev calculation, in which the 
multiple scattering series is summed to all orders. As reference, the dotted line represents twice the two-body total cross 
section 2(j^'^, to which the different second order contributions shown in Fig. [3] are added. The Glauber calculations 
up to second order agrees with the full Faddeev result starting from about 1 GeV. Both, the simplification of the 
Glauber second order term, Eq. (|3.33p . and the high energy limit of Eq. (|2.9p considerably overestimate the exact 
Faddeev result at all energies considered. 



B. Differential Cross Sections 



Similar to the investigations of the total cross section, we want to start comparing the first order calculations of 
the differential cross section at various energies. The differential cross section for elastic n + d scattering is calculated 
from the operator U for clastic scattering, Eq. (|2.3p . and given as function of the laboratory solid angle as 

—— ^ 27r -m 7=== W>dOTo k^dqo • 4.2 

Here xi = cosOi = k'-k represent the laboratory scattering angle for the n + d system, and qo its cm. momentum. For 
calculating the Glauber differential cross section from Eq. (|3.1ip , we have to consider the projection of the momentum 
transfer q on to a plane either perpendicular to the laboratory momentum k or to the sum of incoming and outgoing 
momenta k -I- k', as indicated in Eqs. (|3.24p and (|3.25p . Since both choices are employed in the literature, we will 
consider both in this study. 

Taking the considerations from the previous subsection as guidance, we start our comparison at 500 MeV projectile 
energy and show in Fig. [S] the first order calculations for the differential cross section as function of the laboratory 
scattering angle. The solid line represents the first order Faddeev calculation, and the dashed and dash-dotted lines 
the first order Glauber calculations with the choices (a), Eq. p.24p . and (b), Eq. (|3.25p . for the projection of the 
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momentum transfer. Both choices lead to identical results for small scattering angles. The first small deviation can 
be seen around the first diffraction minimum. Though the Glauber approximation should only be valid in forward 
direction, we plot the entire allowed angle region in order to obtain quantitative insights in the region of angular 
validity. Both Glauber calculations start to deviate from the Faddeev first order calculation as well as from each other 
at 9i K, 50°, which corresponds to a momentum transfer of about 740 MeV. For comparison, we also add the fully 
converged Faddeev result as dash-double-dotted line to the figure. An obvious difference between the first order and 
the fully converged Faddeev calculations is the first diffraction minimum, which is filled in by rescattering corrections. 
Rcscattering contributions are also important for the large angle (high momentum transfer) behavior of the differential 
cross section [T3, [U ■ 

Fig. [S] shows the differential cross section for similar calculations at the considerably higher energy of f.5 GeV. 
Again, the first order Faddeev calculation and the first order Glauber calculations with the two different choices of qj^ 
are close to each other up to 9i « 50°, which now corresponds to a momentum transfer of about f 250 MeV. However, 
the first order calculations differ from the fully converged Faddeev calculation already at much smaller angles, around 
9i K, f2°, which corresponds to roughly 350 MeV in the momentum transfer. While the Glauber formulation is 
intended as high energy approximation, we are interested in exploring how the differential cross sections of the first 
order Glauber and Faddeev calculations compare at lower energies. In Figs. [7] and [5] the differential cross sections for 
laboratory projectile energies of 200 and 100 MeV are displayed. For 200 MeV the two different choices of in the 
Glauber first order calculations already differ at 9i « 25°, corresponding to a momentum transfer of about 270 MeV. 
At roughly the same angle the first order Faddeev calculations also deviates from the converged solution. It is well 
known that the Faddeev multiple scattering series converges very slowly at such low energies [1^, so this finding is 
not surprising. However, at very small angles, all calculations shown roughly agree. This becomes different, if one 
looks at the even lower energy of 100 MeV laboratory projectile energy. Here the first order Glauber calculations 
do not agree with the first order Faddeev calculation even at small angles. It is also very obvious that at this low 
energy the first order Faddeev calculations is quite different from the fully converged one. It is however noteworthy, 
that both first order Glauber calculations are closer to the converged Faddeev result in forward direction than is the 
calculation based on the first order Faddeev term. The two different choices for q^ start to differ around 9i k, 25°, 
which corresponds now only to a momentum transfer of roughly 180 MeV. The overall observation is that at such a 
low energy neither a first order Faddeev calculation nor a Glauber calculation are a good representation of the fully 
converged Faddeev result. 

Next we add the Glauber rescattering correction of Eq. (|3.14p to the differential cross section. In Fig [5] we show 
the effect of this term on the differential cross section at 500 MeV laboratory projectile energy together with the 
contribution of the second order Faddeev rescattering term. While the Faddeev first order rescattering term gives a 
large contribution to the first minimum, the Glauber rescattering term fails to do this, independent of the choice of 
the unit vector fi. Obviously, the minima are more sensitive to the interference of the two terms of the scattering 
amplitude. Thus here the very different structure of the Faddeev and Glauber rescattering terms are clearly visible. 
It should also be noted that in the first minimum one Faddeev rescattering term is already sufficient to coincide with 
the fully converged Faddeev calculation. At slightly larger angles more rescattering terms are required in the Faddeev 
calculation to reach convergence. 

This general feature of the difference between the two approaches does not change with increasing energy. In Fig llOl 
the differential cross section is shown for a laboratory projectile energy of 1 GeV. While the total cross section for 
both, Faddeev and Glauber calculations started to agree at energies of 1 GeV and higher when adding the second 
order correction term in both schemes, this is not the case for the differential cross section. The second order Faddeev 
rescattering term again gives a much larger contribution in the first minimum than the Glauber rescattering correction. 
This trend continues up to 2 GeV, the highest energy we consider here. We also can observe that with increasing 
energy the difference between the two choices of n decreases. To investigate at which energies there are differences, we 
need to consider lower projectile energies. In Fig. [TT] the differential cross section for 200 MeV projectile laboratory 
energy is shown. Here we can first observe, that the second order Glauber correction also increases the cross section 
in the first minimum, but again by far less than the Faddeev rescattering term does. We also observe a distinct 
difference between the two choices of n. At 200 MeV the choice (a) almost has no effect at larger angles, whereas 
choice (b) shows a big contribution. However, since the Glauber expression should be valid only in forward direction, 
one should not put any physical relevance to this difference. 

After having seen that the Glauber rescattering term does not compare well with the corresponding Faddeev term 
in diffraction minima, we want to consider the differential cross section in the very forward direction. This is the 
physical regime for which the Glauber ansatz was developed and where it should perform best. In Fig. [T^] we show the 
three-body differential cross section for very forward angles for laboratory projectile energies of 100 MeV, 200 MeV, 
500 MeV, and 800 MeV. We again compare the results from the first and second order terms in the Faddeev multiple 
scattering series to the fully converged Faddeev calculation, and the first order Glauber calculation with its rescattering 
correction to the Faddeev calculations. Both choices of n coincide for very small angles. Thus we arbitrarily choose 
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(b) for the Glauber calculations of Fig. [121 Since we already established from the comparison of the total cross 
sections, that the Glauber calculations do very well for energies of 1 GeV and higher, we concentrate here on low 
and intermediate range energies. At 800 MeV essentially all calculations agree with each other in the very forward 
direction. At 500 MeV the first order Faddeev term slightly over-predicts the fully converged result. The first Faddeev 
rescattering term has almost no effect, and one needs higher order rescattering contributions to achieve a converged 
result. The Glauber first order calculation differs from the Faddeev first order and only slightly over-predicts the 
converged Faddeev result. The contribution of the Glauber rescattering term is small, and lowers the cross section 
for the very forward angles directly to the full Faddeev result. The situation at 200 MeV is similar. Again, the first 
order Faddeev and Glauber calculations differ, and the Glauber calculation already coincides with the fully converged 
Faddeev result. The Glauber rescattering contribution again slightly lowers the cross section in forward direction. At 
this energy it is interesting to observe that the first rescattering contribution of the Faddeev multiple scattering series 
increases the cross section in forward direction. To reach convergence eight terms of the multiple scattering series 
must be summed up at this energy. At 100 MeV the Glauber first order calculation is already quite close to the fully 
converged Faddeev results, in contrast to the first order Faddeev calculation. The Glauber rescattering correction 
then lowers the cross section essentially to the exact Faddeev result. Thus, for the differential cross section in very 
forward direction the Glauber calculation corrected by the rescattering term captures the exact Faddeev result already 
at 100 McV. Fig. [8] shows that this good agreement is valid to about Oi « 20°. 

The behavior of the differential cross section in the extreme forward direction appears in a very similar fashion in 
the total cross section for elastic scattering, a^i, which is obtained by integrating over the differential cross section. 
Table [IT] lists the total cross sections for elastic scattering as function of the laboratory projectile energy for the 
converged Faddeev calculations, Faddeev calculations including only the first two rescattering terms as well as the 
Glauber calculations. We see that already from 200 McV on the first order Glauber calculations match the converged 
Faddeev results reasonably well. The Glauber rescattering correction is quite small, and always lowers the cross 
section, whereas the first Faddeev rescattering correction increases the cross section, and only the second rescattering 
correction (3rd order in the multiple scattering series) lowers it. It is also obvious from the table that at 100 MeV 
many more terms of the Faddeev multiple scattering series are necessary to obtain a converged total cross section. The 
Glauber result with rescattering correction is similar in size to the Faddeev result with two rescattering corrections, 
though this is most likely accidental. A Glauber calculation for the total cross section for elastic scattering at 100 MeV 
is not trustworthy anymore, since the differential cross sections do not really agree. 

The total cross section is the sum of the total cross sections for elastic scattering and breakup reactions, af^ = 
+ f ■ In Refs. [13, HH this relation was used to estimate the quality of the numerical solutions of the Faddeev 
equation. In a Faddeev calculation the total cross section for breakup reactions is obtained by integration over the 
five-fold differential cross section for breakup, i.e. by summing over all possible breakup configurations. Here we can 
employ this relation to obtain a total cross section for breakup reactions within the Glauber framework. This is done 
in Table Hill The left side lists the Faddeev total cross sections for elastic scattering and breakup reactions with their 
sum, the total cross section, as function of laboratory projectile energy. The right side lists the Glauber total cross 
section obtained from Eq. (|3.12p . The Glauber total cross sections for breakup reactions are obtained by subtracting 
the Glauber total cross sections for clastic scattering from Table HIl from these numbers. For energies less than 500 MeV 
the so obtained breakup cross sections are definitely not competitive with the Faddeev results. However, at 1.5 GcV, 
the Glauber breakup cross section starts to match the Faddeev breakup cross section within 5%. Though this finding 
may look surprising, it could indicate that at sufficiently high energy the dominant ejectiles of a breakup reaction exit 
into the forward cone. That region of phase space however is the one, for which the Glauber ansatz was designed. 
Thus it seems that for GeV projectile energies obtaining a total cross section for breakup reactions from the Glauber 
total and elastic cross section via the optical theorem captures the overall breakup reaction. Of course the description 
of detailed exclusive n + d observablcs will never be accessible in the Glauber formulation. 

V. SUMMARY AND CONCLUSIONS 

In this study we perform fully converged Faddeev calculations for three identical bosons (our model for the n + d 
system) in the intermediate energy regime between 100 MeV and 2 GeV. We calculate total cross sections as well as 
differential cross sections for elastic scattering. The key point of those calculations is the use of vector momenta in the 
formulation, so that all partial waves are automatically included. We then calculated the same cross sections using 
Glauber formulation [25l. l26l. [27| . However a key difference to this early work is that the two-body scattering amplitude 
entering our Glauber expressions is the solution of a two-body Lippmann-Schwinger equation, i.e. the same input as 
is used in the Faddeev calculations. As two-body interaction we employ a superposition of two Yukawa terms, one 
attractive, the other repulsive, for which the parameters are chosen such that a bound state at the empirical deuteron 
binding energy is supported. The deuteron wave function and form factor are also calculated from this interaction. 



10 



By using the identical two-body input for the Faddeev and Glauber calculations we make sure that both calculations 
are based on the same ingredients. This way we can clearly attribute differences in the observables to the different 
treatment of multiple scattering in the two formulations. The first order term in the Glauber formulation is given by 
twice the two-body scattering amplitude, which is folded with the deuteron form factor to give the differential cross 
section and leads to 2(j^^ for the total cross section. The first order Faddeev term is given by T = tP, thus having 
a similar structure. Comparing the total cross sections obtained from the first order term leads to the conclusion 
that from projectile energies of about 200 McV onward both calculations agree. For the differential cross section 
agreement to about 9i w 45° is achieved for projectile energies of 500 McV and higher. For 200 MeV the angular 
range of agreement is already less, and at 100 MeV the first order calculations do not agree with each other. In 
general we can conclude, that in first order the Glauber and Faddeev calculations for observables considered, total 
and differential cross sections for elastic scattering are remarkably close for energies higher than 200 MeV projectile 
energy. 

The consideration of the first rescattering correction in both formulation shows that there is a distinct dependence 
on the observables considered. As a reminder, the Glauber formulation leads only to one rescattering correction, 
whereas the Faddeev formulation as integral equation has an in principle infinite series of rescattering corrections. 
As far as the correction to the total cross section is concerned, the contribution from the first Glauber and Faddeev 
rescattering term starts to become close at 500 MeV and is identical in size at 1 GcV. We also want to point out, 
that this rescattering correction to the total cross section, the shadowing (better named anti-shadowing here, since 
the contribution increases the cross section), is still present at 2 GeV, the highest energy we considered. Thus we 
conclude, that for a good description of the total cross section in the intermediate energy regime the consideration of 
the second order correction in the Glauber as well as in the Faddeev formulation is essential. 

For the differential cross section for elastic scattering the conclusions are twofold. First, when considering the very 
forward direction, we find that at projectile energies of 100 and 200 McV the second order correction brings the 
Glauber calculation of the differential cross section in quite good agreement with the result of the fully converged 
Faddeev calculation. In contrast, the first Faddeev rescattering term even has the opposite effect and moves the 
cross section away from the converged calculation. At projectile energies of 500 MeV and higher, the second order 
corrections in the very forward direction are extremely small and the first order Faddeev and Glauber calculations 
already agree quite well with each other and the converged Faddeev calculation. A similar conclusion can be drawn 
when considering the total cross sections for elastic scattering. Second, when concentrating on the diffractive structure 
of the differential cross section given in our model, we sec distinct differences between the behavior of the Faddeev 
and Glauber second order correction. Whereas the Faddeev correction fills in the minima, the Glauber correction 
can not do that. A similar observation was made in Ref. [l^ in the context of proton elastic scattering from halo 
nuclei considering energies from 100 to 200 MeV per nucleon. However we find that even at much higher energies 
the Glauber second order correction does not come close to the effect of the Faddeev correction. For energies higher 
than 500 MeV the first Faddeev rescattering correction already captures the bulk of rescattering corrections and 
coincides with the converged calculation in the first minimum. Diffraction minima are well known to be sensitive to 
the dynamics of the system or interference effects. The assumption of fixed target particles leading to the specific 
form of the Glauber rescattering correction proves to be too simple to capture the much more involved structure 
of the first Faddeev rescattering term. We also need to point out here, that the diffractive structures shown in the 
differential cross sections of this work are a result of our underlying assumption of identical bosons. The true n + d 
differential cross section has only one minimum and one might speculate that a Glauber calculation at higher energies 
including the rescattering term and being based on nucleon-nucleon (NN) interactions describing the NN observables 
has a chance to describe the n + d total cross and possibly the differential cross section certainly in forward direction 
quite well, provided a relativistic formulation is going to be used. 
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Va [MeV fm] 


MA [fm ^] 


VR[MeV fm] 


Mfl [fm ^] 




626.8932 


1.55 


1438.7228 


3.11 


-2.23 



TABLE I: The parameters and deuteron binding energy for the Malfliet-Tjon type potential (MT3) of our calculation. As 
conversion factor We use units such that ?ic=197.3286 MeV fm = 1. 





Faddecv 


Glauber (a) 


Glauber (b) 


Eiab [MeV] 




[to 


^1 -T^' [fm^] 


o^ei [fm ] 


'^li [fm^] 


[fm^] 


all [fm^] 


1 + 2 rr- 21 

<^eT [fm ] 


100 


40.66 


44.47 


30.27 


26.53 


29.58 


28.24 


32.60 


30.83 


200 


16.81 


18.12 


15.84 


15.31 


15.32 


15.06 


15.74 


15.46 


500 


6.87 


7.04 


6.82 


6.74 


6.59 


6.52 


6.63 


6.56 


800 


4.29 


4.33 


4.34 


4.32 


4.23 


4.19 


4.25 


4.21 


1000 


3.43 


3.46 


3.47 


3.46 


3.44 


3.41 


3.45 


3.43 


1200 


2.91 


2.92 


2.95 


2.96 


2.90 


2.88 


2.91 


2.89 


1500 


2.32 


2.33 


2.41 


2.41 


2.35 


2.33 


2.35 


2.34 


2000 


1.68 


1.68 


1.76 


1.76 


1.78 


1.77 


1.78 


1.77 



TABLE II: The Total cross section for elastic scattering for the first 3 terms of the Faddeev multiple scattering series, the 
fully converged Faddeev calculation, and the first two orders of Glauber calculations based on choices (a) and (b) for the unit 
vector n The superscripts indicate the order in the two-body t matrix. 





Faddeev 


Glauber 


Eiab [MeV] 


[fm^] 


-eT [fm^] 


ar [fm^] 


^Fo'. [fm^] 


aP;(a)[fm2] 


ag.'(b)[fm2] 


100 


34.16 


26.53 


7.63 


29.00 


0.76 




200 


19.00 


15.31 


3.69 


17.39 


2.33 


1.93 


500 


10.30 


6.74 


3.56 


9.49 


2.97 


2.93 


800 


7.22 


4.32 


2.90 


6.77 


2.58 


2.56 


1000 


6.00 


3.46 


2.54 


5.74 


2.33 


2.31 


1200 


5.24 


2.96 


2.27 


5.01 


2.13 


2.12 


1500 


4.37 


2.41 


1.97 


4.22 


1.89 


1.88 


2000 


3.35 


1.76 


1.59 


3.36 


1.56 


1.59 



TABLE III: The total cross section and the total cross section for elastic scattering and breakup reactions as obtained from 
the exact Faddeev calculations are listed in the left part of the table. On the right side the total cross section as obtained from 
the Glauber calculation (including the second order correction) is given. The total cross sections for breakup reactions, ai,r, 
is obtained by subtracting the total cross sections ((a) and (b)) for elastic scattering taken from Table HIl from the total cross 
section a tat- 
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FIG. 1: (Color online) The total cross section (top panel), the elastic total (middle panel) and the breakup total cross section 
(bottom panel) as function of the order t in the multiple scattering series for selected laboratory projectile energies indicated 
in the figure. Starting from the first order in the Faddeev calculation, the next three higher orders (rescattering terms) are 
successively added. 
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FIG. 2: (Color online) The total cross section for three-body scattering as function of the laboratory projectile energy. Shown 
are the first order (in t) Faddeev calculation and twice the two-body total cross section a^^ corresponding to the first order 
term in a Glauber calculation as well as in the expansion of Ref. [2^ . 




FIG. 3: (Color online) The second order correction to the total cross section for three-body scattering as function of the 
laboratory projectile energy. The solid line represents the second order term, tPGotP, of the Faddeev multiple scattering series 
and the long dashed line the correction obtained in Ref. [23 ] from the high energy limit of the 2"'' order Faddeev term. The 
short dashed line corresponds to the 2"'' order Glauber correction Satot from Eq. (|3.32|l . while the dash-dotted line stands for 
the approximation of Eq. p.33p . 
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FIG. 4: (Color online) The total cross section for three-body scattering as function of the laboratory projectile energy. The 
solid line represents fully converged Faddeev calculations, the dotted line twice the two-body total cross section. To this are 
added the correction as extracted from the high energy limit of the 2 order Faddeev term from Ref. [S^l (dash-dotted line), 
the second order Glauber correction 5atot (dash-dash-dotted line) and its approximation (dashed line). 



E, ^ = 500 MeV 

lab 




40 60 



FIG. 5: (Color online) The differential cross section for elastic three-body scattering at 500 MeV laboratory projectile energy 
as function of the laboratory scattering angle. The solid line represent the first order Faddeev calculation. The dashed and 
dash-dotted line stand for the first order Glauber calculations for the two choices (a) and (b) of the direction of the unit vector 
n (see text). For comparison, the dash-double-dotted line represents the full Faddeev calculation. 



16 



R = 1500 MeV 

lab 




ND - 1 




Gl-a - 1 




Gl-b - 1 




ND-full 






40 60 



80 



FIG. 6: (Color online) The differential cross section for elastic three-body scattering at 1500 MeV laboratory projectile energy 
as function of the laboratory scattering angle. The meaning of the curves is the same as in Fig. [5] 




FIG. 7: (Color online) The differential cross section for elastic three-body scattering at 200 MeV laboratory projectile energy 
as function of the laboratory scattering angle. The meaning of the curves is the same as in Fig. [5] 




FIG. 8: (Color online) The differential cross section for elastic three-body scattering at 100 MeV laboratory projectile energy 
as function of the laboratory scattering angle. The meaning of the curves is the same as in Fig. [5] 
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FIG. 9: (Color online) The differential cross section for elastic three-body scattering at 500 MeV laboratory projectile energy 
as function of the laboratory scattering angle. In both panels short dashed line represents the first order Faddeev calculation, 
the long dashed line the first and second order are considered, and the solid line stands for the full Faddeev calculation. The 
left panel shows the first order (dash-dotted line) and first plus second order (double-dash-dotted line) Glauber calculation for 
the choice (a) of the direction of the unit vector h (see text). The right panel contains the same for the choice (b). 
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FIG. 10: (Color online) The differential cross section for elastic three-body scattering at 1000 MeV laboratory projectile energy 
as function of the laboratory scattering angle. The curves have the same meaning as in Fig. (9] 




FIG. 11: (Color online) The differential cross section for elastic three-body scattering at 200 MeV laboratory projectile energy 
as function of the laboratory scattering angle. The curves have the same meaning as in Fig. [9l 
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FIG. 12: (Color online) The differential cross section for elastic three-body scattering in forward direction for 100 MeV (left 
upper panel), 200 MeV (left lower panel), 500 MeV (right upper panel), and 800 MeV (right lower panel) laboratory projectile 
energy. The short dashed line represents the first order Faddeev calculation, the long dashed line the first and second order 
are considered, and the solid line stands for the full Faddeev calculation. The first order and first plus second order Glauber 
calculation for the choice (b) of the direction of the unit vector h are given by the dash-dotted and double-dash-dotted lines. 



